SS-V: 5040 Residual Deformations in an Axially Loaded Plastic Bar

Test No. VNL05 Find residual deformations in a bar axially loaded beyond plasticity.

Definition

Figure 1.

Bar dimensions are 10 x 10 x 200 mm. Strain-stress curve of the bar material is defined by the power law:(1)
σ=Kεn MathType@MTEF@5@5@+=feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaapeGaeq4WdmNaeyypa0Jaam4saiabew7aL9aadaahaaWcbeqaa8qacaWGUbaaaaaa@3C8B@
Where,
K MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKfMBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhiov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaakeaacaWGlbaaaa@3987@
Strength coefficient.
n MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKfMBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhiov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaakeaacaWGlbaaaa@3987@
Must be in the range [0,1].
n MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKfMBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhiov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaakeaacaWGlbaaaa@3987@ =0
Material is perfectly plastic.
n MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKfMBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhiov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaakeaacaWGlbaaaa@3987@ =1
Material is elastic.

The left end of the bar is clamped and the right end is loaded with force F.

The material properties are:
Properties
Value
K MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKfMBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhiov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaakeaacaWGlbaaaa@3987@
530 MPa
n MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKfMBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhiov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaakeaacaWGlbaaaa@3987@
0.26
Poisson's Ratio
0


Figure 2. Corresponded strain-stress curve
Note: Elasticity modulus defined by the first two points of the strain-stress curve is E=2.67324e+10 Pa.

The study was performed for the following load F values: 20000 N, 25000 N, and 30000 N

Reference Solution

One-dimensional analytical reference solution is described here.

At strain ε MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKfMBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhiov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaakeaacqaH1oqzaaa@3A5E@ and stress σ MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKfMBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhiov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaakeaacqaHdpWCaaa@3A7A@ , the residual strain is:(2)
εr=εεe=(σK)1nσE=(F(KA))1n(F(EA)) MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKfMBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhiov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaakeaacqaH1oqzcaWGYbGaeyypa0JaeqyTduMaeyOeI0IaeqyTduMaamyzaiabg2da9maabmaabaWaaSGaaeaacqaHdpWCaeaacaWGlbaaaaGaayjkaiaawMcaamaaliaabaGaaGymaaqaaiaad6gaaaGaeyOeI0YaaSGaaeaacqaHdpWCaeaacaWGfbaaaiabg2da9maabmaabaWaaSGaaeaacaWGgbaabaWaaeWaaeaacaWGlbGaey4fIOIaamyqaaGaayjkaiaawMcaaaaaaiaawIcacaGLPaaadaWccaqaaiaaigdaaeaacaWGUbaaaiabgkHiTmaabmaabaWaaSGaaeaacaWGgbaabaWaaeWaaeaacaWGfbGaey4fIOIaamyqaaGaayjkaiaawMcaaaaaaiaawIcacaGLPaaaaaa@5C95@
Where,
ε MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKfMBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhiov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaakeaacqaH1oqzaaa@3A5E@
Total strain in the bar.
εe MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKfMBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhiov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaakeaacqaH1oqzcaWGLbaaaa@3B48@
Elastic component of the total strain.
εr MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKfMBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhiov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaakeaacqaH1oqzcaWGLbaaaa@3B48@
Residual strain.
A MathType@MTEF@5@5@+=feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKfMBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhiov2DaebbnrfifHhDYfgasaacH8qrps0lbbf9q8WrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaakeaacaWGlbaaaa@3987@
Bar cross-section area.
Then residual displacement at the right end of the bar.(3)
Ur=0L(FK*An   F/(E*A))dx MathType@MTEF@5@5@+=feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaapeGaamyvaiaadkhacqGH9aqpdaGfWbqabSWdaeaapeGaaGimaaWdaeaapeGaamitaaqdpaqaa8qacqGHRiI8aaGccaGGOaWaaOqaa8aabaWdbmaalaaapaqaa8qacaWGgbaapaqaa8qacaWGlbGaaiOkaiaadgeaaaaal8aabaWdbiaad6gaaaGccaGGGcGaaiiOaiabgkHiTiaacckacaWGgbGaai4lamaabmaapaqaa8qacaWGfbGaaiOkaiaadgeaaiaawIcacaGLPaaacaGGPaGaamizaiaadIhaaaa@4EF5@

Results

The bar was modeled as a 3D solid. The left end of the solid was fixed and the right end loaded with axial force (Figure 3).


Figure 3.
The following table summarizes the residual deformations results.
Force F [N] SOL Reference, Residual Displacement [mm] SimSolid, Residual Displacement [mm] % Difference
20000 3.22 3.43 6.78%
25000 9.16 9.172 0.13%
30000 19.077 19.14 0.33%


Figure 4. von Mises stress distribution at load F=30000 N